Swinging and Synchronized Rotations of Red Blood Cells in Simple Shear Flow 
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The dynamics of red blood cells (RBCs) in simple shear flow was studied using a theoretical 
approach based on three variables: a shape parameter, the inclination angle 8, and phase angle <j) of 
the membrane rotation. At high shear rate and low viscosity contrast of internal fluid, RBCs exhibit 
tank-treading motion, where <f) rotates with swinging oscillation of shape and 6. At low shear rate, 
tumbling motion occurs and 6 rotates. In the middle region between these two phases, it is found 
that synchronized rotation of <j> and 8 with integer ratios of the frequencies occurs in addition to 
intermittent rotation. These dynamics are robust to the modification of the potential of the RBC 
shape and membrane rotation. Our results agree well with recent experiments. 

PACS numbers: 87.16.D-, 05.45.-a, 82.40.Bj 
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I. INTRODUCTION 

Soft deformable objects, such as liquid droplets, vesi- 
cles, cells, and synthetic capsules exhibit a complex be- 
havior under flows. For example, in capillary flow,fluid 
vesicles Q, red blood cells (RBCs) @, 11 H S, @,E i, 
and synthetic capsules [9( deform to parachute shapes, 
and RBCs also deform to slipper shapes 0, El IE Q- 
Shape transitions of fluid vesicles occur in simple shear 
flow [l(| EL E3] ■ Membrane wrinkling appears for fluid 
vesicles after inversion of an elongational How Il3j, Il4j 
and for synthetic capsules in simple shear flow [la . 1 1 €31 ] . 
Among these soft objects, RBC have received a great 
deal of attention, since they are important for both fun- 
damental research and medical applications. In microcir- 
culation, the deformation of RBCs reduces the How resis- 
tance of microvessels. In patients with diseases such as 
diabetes mellitus and sickle cell anemia, the RBCs have a 
reduced deformability and often block the microvascular 

flow @, [IE [IE Hill. 

In a simple shear flow with flow velocity v = jyex, 
fluid vesicles and RBCs show a transition from a tank- 
treading (TT) mode with a constant inclination angle 9 
to a tumbling (TB) mode with increasing viscosity of the 
internal fluid rji n 0, [22|, [IE [24[ or membrane viscosity 
Vmb [IE EH • This transition is described well by the the- 
ory of Keller and Skalak (KS) [2l[, which assumes a fixed 
ellipsoidal vesicle shape. Experimentally, synthetic cap- 
sules and RBCs show the oscillation of their lengths and 
0, called swinging (SW) [H, HE EH , during TT motion, 
and RBCs also transit from TB to TT with increasing 7 
[lEH?]]- Recently, this dynamics was explained by the KS 
theory with the addition of an energy barrier for the TT 
rotation caused by the membrane shear elasticity [26l. l28l] . 
More recently, this transition was also obtained by simu- 
lations [IE EH- However, the detailed dynamics has not 
yet been investigated. 
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For fluid vesicles in high shear flow, shape transitions 
fToL [Til . Il2| occur, and a swinging phase O, EiE EiE HE 
|34| . where the shape and oscillate around ~ 0, ap- 
pears between the TT and TB phases. This SW mode is 
also called trembling 0, EiE ESj] or vacillating-breathing 
f32| : it is explained by the KS theory extended to a de- 
formable ellipsoidal vesicle [3l| and the perturbation the- 
ory for a quasi-spherical vesicle [H, [H, EH • Shape defor- 
mation plays an essential role in the SW of fluid vesicles. 
The deformation is not necessary to explain the SW of 
elastic capsules [H, HE HI HE Hi HE ES EH but is re- 
quired for quantitative analysis. In this letter, we extend 
the theory in Ref. [28| to include the shape deformation of 
RBCs and investigate the dynamics of deformable RBCs. 

The internal fluid of RBCs behaves as a Newtonian 
fluid since RBCs do not have a nucleus and other in- 
tracellular organelles. The RBC membrane consists of 
a lipid bilayer with an attached spectrin network as cy- 
toskeleton. The lipid bilayer is an area-incompressible 
fluid membrane. The shear elasticity of the composite 
membrane is induced by the spectrin network. Under 
physiological conditions, an RBC has a constant volume 
V = 94^m 3 , surface area S — 135^m 2 , rj m = 0.01Pa-s, 
?7mb ~ 1CP 7 — 10~ 6 Ns/m, membrane shear elasticity 
/i = 6x 10~ 6 N/m, and bending rigidity k = 2x 10 _19 J 

dEl, El mil. 

The models and results arc presented with dimension- 
less quantities (denoted by a superscript *). The lengths 
and energies are normalized by Rq = yj S/4tv and (j,Rq, 
respectively. For RBCs, they are Rq = 3.3 /im and 
fiRl = 6.5 x KT 17 J. There are two intrinsic time units: 
the shape relaxation time r = 770-R0/M by the shear elas- 
ticity n, and the time of shear flow I/7; the reduced 
shear rate is defined as 7* = jr. The relative viscosi- 
ties are 77^ = i]- m /r]o and ?7* lb = i] mh /r] R( h where rj is 
the viscosity of the outside fluid. In typical experimental 
conditions, the Reynolds number is low, Re< 1; hence, 
the effects of the inertia are neglected. 

In Sec. [Ill we describe the extended KS theory [28[ for 
an elastic capsule with a fixed ellipsoidal shape, and the 
phase behavior of the capsule. In Sec. IIIII we introduce 
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the shape equation for deformable RBCs and present the 
dynamics of deformed RBCs. The dependence of the 
function shape of the RBC free-energy potential is de- 
scribed in Sec. HVl Discussion and summary are given in 
Sec. [V] and Sec. IVI1 respectively. The comparison with 
experimental results is presented in Sec. |V] 



II. DYNAMICS OF ELASTIC CAPSULES WITH 
FIXED SHAPE 
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Models 



1. Keller-Skalak Theory 



Keller and Skalak (KS) [2l| analytically derived the 
equation of the motion of vesicles or capsules based on 
Jeffery's theory [39[. In the KS theory, the vesicles are 
assumed to have a fixed ellipsoidal shape, 
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where a,\ are the semi-axes of the ellipsoid, and the coor- 
dinate axes Xi point along its principal directions. The 
Xi and X2 axes, with a± > a 2 , are on the vorticity (xy) 
plane, and the X3 axis is in the vorticity (z) direction. 
The maximum lengths in three directions are L\ — 2a\ 1 
L 2 = 2a 2 , and L3 = 203. The velocity field on the mem- 
brane is assumed to be 
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The energy W ex supplied from the external fluid has to be 
balanced with the energy dissipated in the vesicle, W ex = 
-Din+-Dmb7 where D m and -D m b are the energies dissipated 
inside the vesicle and on the membrane, respectively. The 
motion of the vesicle is derived from this energy balance. 
Then the motion of the inclination angle 9 is given by 
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The membrane- viscosity term has been derived by Tran- 
Son-Tay et ai 17]. The factors appearing in Eqs. ([3][5]) 
are given by 

fa = 2/(01/02 + 02/01), 
/1 = 0.5(oi/a 2 - a 2 /ai), 



where E s is an integral over the membrane surface, and 
n is the normal vector of the surface. 

For B > 1, a stable fixed point 9 = 0.5 arccos(l/i?) 
exists, and TT motion occurs, while for B < 1, there is no 
fixed point, and the angle 9 periodically rotates (TB). As 
77i* n or ?7j^ b increases, the transition from TT to TB motion 
occurs, where B decreases from B > 1 to B < 1. The 
membrane viscosity ?7 m b and the internal viscosity f]* n 
have a similar effect; hence, an effective internal viscosity 
can be defined as ri* s — r]* n + fzri*^ h . The factor / 3 in 
r]* s depends on the vesicle shape and can give different 
dynamics for deformable vesicles, in particular for shape 
transformations between prolate and oblate vesicles [10, 

The KS theory quantitatively predicts the TT-TB 
transition with increasing -q* s . However, it cannot ex- 
plain the TB-TT transition with increasing 7. In the KS 
theory, vesicle motion does not depend on 7 except that 
the TT or TB rotation velocity increases linearly with 7. 



2. KS Theory with an Energy Barrier 

Skotheim and Secomb extended the KS theory to take 
into account an energy barrier during TT membrane 
rotation [28|. For RBCs and synthetic capsules with 
non-spherical rest shape, their membranes are locally 
deformed during the TT rotation. Fischer experimen- 
tally demonstrated that the RBC membrane rotates back 
to the original position when the shear flow is switched 
off [!(|. To describe the energy barrier, a phase angle <ft 
and free energy potential F{4>) are introduced; see inset 
of Fig[TJ The potential is periodic, F{<p + nn) = F(<j>) 
and <p = at the rest shape. Thus, the motions of the 
inclination angle 9 and phase angle 4> are given by 
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where Co = 3/ 2 /87r/i. The equations of the original KS 
theory are recovered in the absence of barriers of the 
free energy F, i.e., dF*/d(j> = 0, where u — d<j>/dt is 
independent of (f>. 



FIG. 1: (Color online) Dynamic phase diagram of elastic 
capsules with fixed shapes calculated from Eqs. (|6]) and 
with the potential Fq = Eo sin 2 ((f)). (a) Viscosity r)*^ depen- 
dence for the RBC-like oblate shape with L2/L1 = 0.25 and 
L3/L1 = 1. (b) Aspect ratio L2/L1 dependence at the TT- 
TB transition viscosity ry* ff = r]l of the KS theory (U e = 0). 
The boundary lines of the TB (brown), TT (violet), and syn- 
chronization regions [/ r 9 ot : ff ot = 1:1 (red), 2 : 1 (blue), and 
3 : 1 (green)] are shown. 



Skotheim and Secomb used a simple potential Fq(<P) = 
EQS\v?((j)) and a reduced energy U c — /2-E0/ 2/1777^ = 
Eq c oH*V* ■ We employ the potential F in this section 
and describe the dependence on the potential shape in 
Sec. [TV] Eqs. ([6]) and (O are numerically integrated us- 
ing the fourth-order Runge-Kutta method. An oblate 
capsule with L 2 /Li — 0.25 and L 3 /Li = 1 is used as a 
model RBC. 



B. Results 

Figures [10 show the dynamics of elastic capsules with 
fixed shapes. At low shear rate 7 (large U e ), the capsules 
show TB motion, since the free energy barrier locks the 
phase angle at (j> ~ 0. In TB, 9 rotates but <\> osculates; 
see the trajectory at U c = 0.81 in Fig. [3] At higher 7 
(smaller i7 e ) and low rji n , TT motion occurs. In TT, 9 
oscillates (swing) , and <fi rotates instead of 9; see the tra- 
jectory at U c = 0.76 in Fig.[3J The oscillation of (ft or 9 in 
TB or TT occurs with the rotation frequency f? ot or ff ot , 
respectively. Here, an angle change of 7r is counted as one 
rotation. Skotheim and Secomb [28j reported an inter- 
mittent phase between the TT and TB phases, where 
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FIG. 2: (Color online) Rotation frequency / rot of the incli- 
nation angle 9 and phase angle <f> for the fixed oblate shape 
with L2/L1 = 0.25 at ?7* ff = 1. Numbers represent / r e ot : ff ot . 



both rotations of 9 and cj) occur. The 4> (TT) rotation is 
intermittently interrupted by the 9 (TB) rotation slightly 
above the maximum energy barrier J7* of the TT phase. 
However, we found that the phases of synchronized rota- 
tion of 9 and <f> also exist in this middle range between the 
TB and TT phases. An infinite number of synchroniza- 
tion phases with integer ratios of / r 8 ot and ff ot exist; see 
Fig. O This type of synchronization is called the Devil's 
staircase [4lj |. 

The trajectories of the synchronized rotations with 
/rot : /rt,t = 1 : 1 and 2 : 1 are shown in Fig. G3 The 
former has the widest U c range. The approach to syn- 
chronization is explained by the return map 4> n +i{4>n) 
at 9 = it/2 — nir in Fig. 2J The curve n +i(0„) shifts 
to the left with increasing U e and has two (stable and 
unstable) crossing points with the line <p n +i = 4>n at 
0.77823 < U c < 0.78722, where the capsules approach 
the limit cycle with 4> n +i — 4> n . At slightly below or 
above the synchronization region (0 < 0.77823 -C/ C <1 
or < U c — 0.78722 <C 1), intermittent rotation appears, 
like near the TT and TB regions [28[; see the trajectory 
at U c = 0.7782 in Fig. H 

The qualitative behavior of capsules does not depend 
on the aspect ratios L 2 /Li and L 3 /Li. However, the 
ranges of the synchronized rotations are narrower for 
more spherical capsules, as shown in Fig. [IJb). Thus, 
it would be difficult to observe synchronized rotations 
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FIG. 3: (Color online) Dynamics of oblate capsules for the 
fixed shape. Top panel: time development of angles 9 and cf> 
at U c = 0.78. Middle and bottom panels: trajectories on the 
phase space (8, <j>). The capsules exhibit TT and TB rotations 
at U c = 0.76 and 0.81, respectively. Synchronized rotations 
with / r s ot : ff ot =1:1 and 2:1, and intermittent rotation are 
observed at U e = 0.78, 0.793, and 0.7782, respectively. The 
other parameters are the same as in Fig. [3] 




FIG. 5: (Color online) Effective force <9F*(«i3, 4>)/d(f> at <f> — 
and 7r/2. Solid and dashed lines represent the simulation 
data and fit functions, respectively. In the inset, a snapshot 
of the model RBC elongated by mechanical forces is shown. 

in quasi-spherical capsules. Note that the dynamics is 
independent of L3/L1 at rj* s = r)* (B = 1). 

A fluid vesicle has no membrane shear elasticity (U c = 
0), and shows no synchronization between 9 and <j) since 
the potential does not depend on <fi; see Fig. [TJa). At 
< ?7* ff — 77* <C 1, intermittent 9 rotation also appears 
for the fluid vesicle. The tumbling frequency is give n by 
/£ t = 7VI - B 2 /2tt (x - T]* in the KS theory [Hf . 

A similar frequency dependence on U c is obtained in the 
extended KS theory: / r e ot oc \JU C — slightly above 

III. DYNAMICS OF DEFORMABLE RBCS 



A. Shape equation 




FIG. 4: (Color online) Return map of oblate capsules for 
the fixed shape. The angle <f> at 6 = 7r/2 — mr is plotted. 
Synchronized rotation with / r e ot : ff ot = 1:1 occurs at 
0.77823 < (7 C < 0.78722. A dashed line represents <f) n +i = <f> n - 
The other parameters are the same as in Fig. [2] 



Previously, we extended the KS theory to include the 
shape deformation of fluid vesicles, on the basis of the 
perturbation theory [H, [H, H3] °^ quasi-spherical vesicles 
[3ll | . It showed very good agreement with experimental 
data [Hj]. Here, we have adapted it to RBC dynamics. 
The shape parameter a i2 = {Li — L 2 ) / {L\+ L 2 ) does not 
increase monotonically for elongation, because of RBC 
dimples. Therefore, the shape parameter 0:13 = {L\ — 
L^)/{Li + L3) is employed, where da\3 / dai2 = 2 for an 
oblate ellipsoid (0:13 — 0). The equation of the shape 
evolution is given by 
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where A = 45/2tt(32 + 23^ + 16r]^ h )V* and A x = 
60/(32+23r?i* n +16?4 b ). Here, the terms of tt^k are added 
in Aq and A\ based on the theory in Refs. [33j, [34{. This 
revision improves the r)* lh dependence of fluid vesicles in 
Ref. [31(; the phase diagram for r?j^ b becomes similar to 
that for r]* n . 

The free energy F(ai3,0) is estimated by the sim- 
ulation of a model RBC with <f> — 0,7r/2 elongated 



5 



0.04 



0.03 



0.02 - 



0.01 



I I I I I I I I 

I (a) 

- Tank-Treading 0:1 / 


/ / ' 

7 1:1/ 


ii i i ii i 

Ay // 

3:1 " 


- ^^.^^^^ Tumbling 1 :0 
i i i i I i i i i I i i i i 



0.5 ri* 1 

' in 



1.5 



0.05 

^* 

o 



*\ 1 


i i i 


1 1 " 

(b) - 








, , i 


i i i 


i A 



0.016 



0.017 



Y 



0.018 



FIG. 6: (Color online) Dynamics of the deformable RBCs in 
the simple shear flow calculated from Eqs. ([6]), (0, and (JSj . 
(a) Dynamic phase diagram, (b) Rotation frequency / ro t of 9 
and (f> at rj* n = 0.5 {<q = 0.02Pa-s ). 



by mechanical forces, where the RBC membrane is 
modeled as a triangular network [12 ]. In the simula- 
tion, 578 vertices are connected by a bond potential 
t/bond = (&i/2)(r - r ) 2 {l + (k 2 /2)(r/r - l) 2 } with 
fj, = (V3/4)fci = 6 x 10- 6 N/m, k=2x 10~ 19 J, and 
k 2 = 1. The area and volume of the RBC are kept con- 
stant by harmonic potentials. Our simulation reproduces 
the force- length curves of the optical- tweezers experiment 
[43| and other simulations [H, 0, Hf|; see Fig. 5 in 
Ref. [l2[. The effective force dF{a.xz,4>)/d(j) is estimated 
from these force-length curves; see Fig. O The model 
RBC at 4> — 7r /2 has 9 x 10~ 18 J higher energy than at 
= with 013 = —0.1 in the absence of external forces. 
This height of the e nerg y barrier agrees with the value 
E Q = 10- 17 J in Ref. Jgi], which was estimated from Fis- 
cher's experiments [40J. Abkarian et al. estimated the 
height as E = fj>S(a 1 /a 2 - a 2 /a 1 ) 2 /2 - 3 to 7 x lO" 10 
[m 2 ] x/i based on the velocity field of the KS theory 
[Eq. @] [26]. However, it gives much higher barrier 
Eq ~ 10~ 15 J for (j, — 6 x 10~ 6 N/m, or smaller shear 
modules n ~ 10- 8 iV/m for E a = 10~ 17 J. Since the KS 
velocity field does not satisfy the local area conservation 
of membrane, it may give more stress on the membrane 
than the area-conserving velocity field [46| , and the bar- 
rier height may be overestimated. 

In this section, we employ the free-energy potential 




FIG. 7: (Color online) Shear rate 7* dependence of the de- 
formable RBCs at rj* a = 0.5. (a) Rotation frequency / ro t of 9 
and (j>. The maximum and minimum of (b) the shape param- 
eter 0:13 , (c) 9, and (d) cj> are shown. 
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Fi (0:13) + ^2(0:13) sin ((j)). The dependence 



on the potential function is discussed in Sec. IIVI Instead 
of an interpolation [Hi EH [HH , we used fit functions to 
obtain smooth functions for the numerical calculations: 
The normalized potentials F*(ai 3 ) = 5a 2 3 + (40/ 3)af 3 + 
(230/4)af 3 , (cm) = 0.2 + 0. 8013; the shape parameter 
a 12 = 0.56 + 0.35ai3-0.23a 2 3 + 0.034a? 3 ; the coefficients 
f 2 = 0.6018 + 0.064cv 13 -0.19a 2 3 - 0.42af 3 , f 3 = 0.734 + 
0.54ai3 + 0.91a 2 3 + 3.2af 3 ; and off* = 0.7. Eqs. ©, 0, 
and ([5]) are numerically integrated using the fourth-order 
Runge-Kutta method. 

In this model, the viscosity ratio of the membrane and 
inner fluid and the reduced volume are fixed at ^mb/^fn = 
3.1 and V* = V/(iirR%/3) = 0.64. Experimentally, the 
viscosity 770 of outside fluid is typically varied and r] m h 
and ?7i n are fixed at physiological values. Thus, r/* n = 1 
corresponds to 770 = 0.01Pa-s. The reduced shear rate 
7* = 1 corresponds to 7 = 90s -1 or 180s -1 , at rj* a = 0.5 
or 1, respectively. 



B. Results 

The phase diagram and rotation frequencies of the 
deformable RBC in simple flow are shown in Fig. [5] 
The shape deformation does not qualitatively change the 
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FIG. 8: (Color online) Dynamics of the deformable RBCs at 
rj* n — 0.5. (a) Time evolution of ai3, 0, and <f> at 7* = 0.05. 
(b) Trajectories in the phase space (Q13, 0) at 7* = 0.02, 
0.05, and 0.1. (c) Shear rate 7* dependence of the phase 
difference between Q13 and 9 oscillations at rj* n = 0.1. Solid 
and dashed lines represent ^ ax /Ttt and tfj\ lxx /Tu, respectively, 
where T tt is the TT rotation period, and £^ ax (tA in ) are the 
time difference between the maximum (minimum) of 013 and 
9; see the explanation in (a). 



phase diagram. The tumbling phase appears in the bot- 
tom of Fig. EJa), since 7* is used for the vertical axis 
instead of U c oc 1/7*. 

In the TT phase, the shape parameter a.13 and 8 os- 
cillate with the frequency }f ot ; see Figs. [TJb) and[8ja). 
As 7* increases, the time-average {a\z) increases, and 
the SW oscillation amplitudes of 0113 and 8 decrease; see 
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FIG. 9: (Color online) Time averages and peak-to-peak 
amplitudes of the (a) 0:13 and (b), (c) 9 oscillations during 
TT motion. Solid lines represent the deformable RBCs with 
the potential To at r/* n — 0.1, 0.5, and 0.8. Dashed-dotted 
lines represent the deformable RBCs with the potential T a at 
r/* n = 0.5. Dashed lines in (b) and (c) represent the RBCs 
with the fixed shape at »7* ff = 0.5 and 2. The shear rate 7 is 
normalized by the minimum shear rate jtt of the TT phase. 



Fig. 02a) . The peak-to-peak amplitude # a mp is inversely 
proportional to 7* for both the deformable and fixed- 
shape RBCs. For the fixed-shape RBCs, the mean angle 
(8) increases with increasing 7* and reaches the angle of 
the KS theory at U c = (7* — * 00). For the deformable 
RBCs, (8) has a maximum and then decreases because 
of the elongation of the RBCs; see Figs. E£b). 

In the TT phase, ai3 and 8 oscillate with a fixed phase 
difference; see Fig. [5] The phase difference is calculated 
from the time difference i^ ax between the maximum val- 
ues of CV13 and 8 and from i™ n between the minimum 
values. The difference between i^ ax /T t t and t^ in /T tt rep- 
resents the asymmetry of the oscillation functions. In a 
sinusoidal function, i^ ax = tf£ m . It is found that a±3 and 
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FIG. 10: (Color online) Time averages and peak-to-peak 
amplitudes of the (a) a\z and (b) (f> oscillations during TB 
motion. Solid lines represent the deformable RBCs with the 
potential Fo at -q* n — 0.1, 0.5, and 0.8. Dashed-dotted lines 
represent the deformable RBCs with the potential F a at ri* n = 
0.5. Dashed lines in (b) represent the RBCs with the fixed 
shape at r/* s — 0.5 and 2. The shear rate 7 is normalized by 
the maximum shear rate 7tb of the TB phase. 



The angle ip is the rotational angle ip = arctan(cCj/y,), 
where r; = (xi,yi,Zi) is the position of a tracer on 
the membrane. The phase angle <p is defined as (j) — 
aictainixiLz/yiLi). The angle ip can be defined without 
assuming an ellipsoidal shape. The potential F a has the 
sharpest peak. 

As the peak of potentials (Fo, F4, F a ) sharpens, the 
transition shear rates 7 t * b and 7^. increase (U c decreases), 
since the maximum forces dF/d(j) increases at constant 
Eq. Two limit cycles of the synchronized rotations can 
coexist for F4 and F a , while no coexistence is observed 
for Fq. The capsules approach different limit cycles with 
increasing and decreasing 7*. For the synchronized ro- 
tations with / r e ot : ff ot = 1 : 1, a steeper but continuous 
change appears in the / ro t curve for F4 than for Fq, and a 
discrete change with hysteresis appears for F a ; compare 
Figs. d^a) and Figs. Efb), (c). 

The coexistence of two synchronized rotations oc- 
curs for the deformable RBCs with the free energy 
F a (ai3» = Fi(ai 3 ) + F 2 (a 13 )sm 2 (<p). The shape pa- 
rameter ai3 and 9 steeply change at <f> ~ ±0.57r during 
TT motion, as shown in Fig. [T2j The oscillations of «i3 
and 9 do not approach sinusoidal curves, unlike the case 
for F - 7 compare Figs.^c) andd^b). In the TT and TB 
phases, the averages and amplitudes of the angles 9 and 
4> do not show significant differences between the poten- 
tial shapes in <f> {Fq, F4, and Fa), and in ai3 (fixed or 
deformable shape); see Figs. I§1 and [TU1 



9 show in-phase oscillation {t&/T t t — 0) at small 7*, and 
the phase difference approaches 7r/4 with increasing 7*. 
Walter et al. experimentally observed a phase difference 
of 7r/4 for synthetic capsules; see Fig. 7 in Ref. [l5j ]. Our 
results agree with their experiments. 

In the TB phase, (0) decreases and </> amp increases with 
increasing 7* ; see Fig. [TDJ When the energy barrier at 
<f)/n = —0.5 is overcome at 7 t * b , <j> begins to rotate. The 
average (aia) decreases, since F (ai3,4>) nas a minimum 
at ai3 < for (f> =/= 0. In TB, there is no significant 
difference between the deformable and fixed-shape RBCs. 



IV. DEPENDENCE ON POTENTIAL 
FUNCTION 



We compared the dynamics of RBCs with the fixed 
and deformable shapes in the previous section. The fix- 
ation of the RBC shape can be interpreted as the bend- 
ing rigidity k — > 00, or F(ai3,<f) = kai3 2 + F<t>{4>) with 
k — ► 00. This difference of the potential functions in 
ai3 does not change the dynamics greatly. In this sec- 
tion, we investigate the dependence on the potential func- 
tions in <fi. For fixed-shape capsules, we compare the dy- 
namics with three potential functions: Fq — Fosin 2 (0), 
F 4 = F o sin 4 (0), and F a = F sin 2 ^); see Fig. ITTTa). 



V. DISCUSSION 



Let us compare our results with the experimental re- 
sults The TT-TB transition is observed at 7* = 
0.01 - 0.04 for 7]* n = 0.45 in Ref. [26]. It is in good 
agreement with our results, 7 t * b = 0.0158455 and 7 t * t = 
0.0173455 with F and 7 * b = 0.0272 and 7* = 0.028 
with F a . Abkarian et al. also reported the dependences 
in the TT phase, / ro t °c 7 (T tt oc 1/7) and the SW am- 
plitude # amp cx I/7. We also obtained these dependences 
for all types of potentials; see Figs.[7£a,) and[5Jc). Most 
of the dynamics in our calculations are qualitatively in- 
dependent of the potential shapes, and the quantitative 
differences are smaller than the distribution widths of 
the experimental data. In the experiments, the transi- 
tion shear rate and the other quantities have wide dis- 
tributions because of the polydispersity of RBCs. RBCs 
become smaller and more viscous with age [l7l [l8j . The 
viscoelasticity of the RBC membrane is changed by some 
diseases such as diabetes mellitus [lj|. Further experi- 
mental and simulation studies are needed to tune up the 
RBC potential in our theoretical model. 

In experiments, it is difficult to distinguish the inter- 
mittent rotations from transient rotations [26[. The in- 
termittency has not yet been obtained by numerical simu- 
lations [29|,[3(|. The synchronized rotations, in particular 
with / r e ot : ff ot = 1:1, would be much easier to observe 
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FIG. 11: (Color online) (a) Three potentials Fq = Eo sin 2 (0), 
Fi — Eo sin 4 (</>), and _F a = _Eo sin 2 (<^) . The rotation frequen- 
cies / r ot of 8 and for fixed-shape capsules with (b) F4 and 
(c) F a at ?7* ff = 1. Dashed or solid lines in (b) represent data 
obtained with increasing or decreasing 7*, respectively. The 
other parameters are the same as in Fig. [2] 



in experiments and simulations. Kessler et al. [29( ar- 
gued that the intermittent phase might be an artifact of 
the theory, since they did not obtain it in their simula- 
tions. In our study, however, the modifications of the 
theory do not qualitatively change the phase diagram. 
The Devil's staircase is the general dynamics on a torus 
(two-dimensional plane with periodic boundary condi- 
tion). We believe that the intermittent and synchronized 
rotations occur at higher 7 than their simulated values. 

Interestingly, for TT motion slightly above 7 t * t , the an- 
gle 9 oscillates crossing 6 — 0; see Figs. [7{c). Fluid vesi- 
cles also show SW oscillation of the vesicle shape and 9 
around 9 ~ between the TT and TB phases by a dif- 
ferent mechanism, where the shrinkage of the vesicles at 
9 < induces a change of the dynamic mode from TB 
(B < 1) to TT (B > 1) in the generalized KS theory [Uj]. 
Previously, we distinguished that in the SW oscillation 
of the elastic capsules, 9 is always positive and the shape 
deformation is negligibly small, while in SW of fluid vesi- 
cles, 9 changes its sign and the shape shows large defor- 
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FIG. 12: (Color online) Dynamics of the deformable RBCs 
with the potential _F a . (a) Time evolution of 0:13, 8, and <j> at 
7* = 1 and 77i* n = 0.5. (b) Solid and dashed lines represent 
tTVTtt and tf n /T tt , respectively. 



mation. However, we know now that the condition for 9 
is not always true. The clear difference is the dependence 
on 77j* n or 7j^ b . SW induced by the shape deformation ap- 
pears only in a narrow range of the viscosity, whereas 
SW induced by the membrane shear elasticity appears 
at a wide range of the viscosities with no lower viscosity 
limit. In the future, it will be interesting to investigate 
the coupling of different oscillation mechanisms in elastic 
capsules. 



VI. SUMMARY 



In summary, we described the dynamics of RBCs in 
simple shear flow using a simple theory. The phase 
diagram of RBCs is divided into three regions: tank- 
treading, tumbling, and intermediate regions. In the in- 
termediate regions, RBCs exhibit intermittent or syn- 
chronized rotations of the inclination angle 9 and phase 
angle 0. Synchronized rotations, in particular with / r e ot : 
ff ot = 1:1, would be much easier to experimentally ob- 
serve than intermittent rotations. In the TT (TB) phase, 
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the shape and 9 (4>) oscillate with the frequency of <f> (9) 
rotation. The coexistence of two synchronized rotations 
can appear when the potential function of <f> has a sharp 
peak. The other dynamic properties are not sensitive 
to the function shape of the free-energy potential. We 
focused on the dynamics of RBCs in this paper, but the 
resulting dynamics would be generally applicable to other 
clastic capsules. 
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